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Solutions of Ginzburg-Landau eqns. coupled with three 
dimensional Maxwell eqns. reveal intriguing magnetic re- 
sponse of small superconducting particles, qualitatively differ- 
ent from the two dimensional approximation but in agreement 
with recent experiments. Depending on the radius and thick- 
ness first or second order transitions are found for the normal 
to superconducting state. For a sufficient large radius of the 
disc several transitions in the superconducting phase are ob- 
tained which correspond to different angular momentum giant 
vortex states. The incorporation of the finite thickness in the 
calculation is crucial in order to obtain agreement with the 
position and the size of these jumps, and the line shape and 
magnitude of the magnetization curves. 
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Recent advances in micro-fabrication technology and 
measurement techniques have allowed the first studies of 
thermodynamic properties of well controlled mesoscopic 
superconducting particles The samples are meso- 

scopic in the sense that their size is comparable to the 
Ginzburg-Landau (GL) coherence length. Buisson et al 
performed magnetization measurements on an ensem- 
ble of discs with large separation between them in order 
to make the dipolar interaction between the discs negli- 
gible. They found oscillatory behavior in the magneti- 
zation near the superconducting transition temperature 
and showed that the linearized Ginzburg-Landau (LGL) 
eqn. are able to explain qualitatively part of their ex- 
periments but there are some major discrepancies in size 
and position of the jumps in the magnetization. Recently 
Geim et al |3| used sub-micron Hall probes to detect the 
magnetization of single superconducting discs with size 
down to 0.1 jjLm. At different applied fields the discs 
show various kinds of phase transitions within the super- 
conducting state and between the superconducting and 
normal state which can be first or second order depend- 
ing on the sample dimensions and temperature. The aim 
of this letter is to explain this intriguing behavior and 
to give a quantitative analysis of these magnetization ex- 
periments. 

A number of earlier works studied: 1) disc geometries 
in the framework of the LGL eqn. where a uniform mag- 
netic field is assumed [|l|,||j4[|^ in the discs, or 2) cylin- 
drical geometries in the framework of the non-linear GL 
eqns. P|. The first type of approximations is reasonable 
if one is interested in the superconducting-normal state 



boundary where the order parameter is very small and 
the magnetic field equals the external one. The second 
type of approach is also valid away from this boundary 
but does not work for discs of finite thickness. Here we 
are interested to explain the full magnetization curve as 
a function of the external magnetic field for mesoscopic 
discs with finite thickness. Magnetization is the ideal tool 
to understand the superconducting state deep inside the 
phase boundary. 

It is known that a type H superconducting cylinder in a 
magnetic field parallel to its axis can exist in three phases 
of superconductivity. Below Hd we have a pure super- 
conducting state, between Hd and Hc2 there is the mixed 
state and between Hc2 and we have surface super- 
conductivity or the giant vortex state 0. As the height 
of the cylinder is reduced, so that it becomes a disc, the 
magnetic field penetration into the sample is determined 
by the penetration length A as well as the disc thickness, 
due to geometrical form factors. This makes the above 
simple divisions no longer applicable. Our problem re- 
quires a 3D solution, instead of their 2D version which 
turns out to be essential in order to understand the ex- 
periments of Refs. and 

We consider a superconducting Al discs (k=.28) with 
radius R and thickness d immersed in an insulating 
medium. Thin film discs are known to behave as type 
II samples |^ and can be described by the GL theory 
For mesoscopic Al samples (squares and thin wires), 
the GL theory has been successfully employed to explain 
the phase boundary Hence as a first approximation, 
neglecting the non-local effects, we solve the system of 
two coupled GL eqns. 
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where the density of superconducting current j is given 
by 



J= f5'*V*-*V**) - — |*|^A (3) 
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The boundary conditions on the disc corresponding to 
zero current density in the insulator medium is 
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where the subscript n denotes the component normal to 
the disc surface. The boundary condition for the vector 
potential is such that far away from the superconduct- 
ing disc the field equals the appHed field H = (0, 0, Hq), 
i.e., A\^^ao = Here denotes the azimuthal 

direction and p is the radial distance from the disc center. 

Using dimensionless variables and the London gauge 
divA = 0, we rewrite the system of eqns. (0^) into the 
following form 
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Here the distance is measured in units of the coher- 
ence length ^ = h/y/— 21710, the order parameter in 
■00 = \/— a/ f3, the vector potential in c?i/2e^, k = A/^ 
is the GL parameter, and A — Cy/rnpK / Aeipo is the 
penetration length. We measure the m agnetic fi eld in 
Hc2 = ch/2eS,'^ = kV2Hc, where He = i/-47ra//9 is the 
critical field. The difference of the Gibbs free energy G 
between the superconducting and the normal state mea- 
sured in H'^V/8tt can be expressed through the integral 
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over the disc volume V = irR^d, where Aq — e^Hop/2 
is the external vector potential, and j — (5'*V5' — 
^'V^'*)/2i — I^PpA is the dimensionless superconducting 
current. 

For thin discs, the magnetic field is uniformly dis- 
tributed along the z— direction. When the disc thick- 
ness becomes comparable to the penetration length, the 
magnetic field is expelled from the disc due to the Meiss- 
ner effect. The field penetrates only a depth A inside 
the disc. Therefore, the variation of the vector po- 
tential in the direction parallel to the applied field be- 
comes rather strong for d > X. Nevertheless, this does 
not lead to essential variations of the order parameter 
in this direction, in discs thinner than the coherence 
length. Representing the order parameter as a series 
= '^f^cos{k'Kz/d)'^k{p), which obeys the bound- 
ary condition (jj) at the disc sides z — ±c?/2 and using 
the first GL eqn. one can verify that the part of 

the order parameter which is uniform in the z-direction 
i.e., ^'o gives the main contribution to the expansion for 
(TT^/d)^ ^ 1. Therefore, we may assume a uniform order 
parameter along the direction and average the first 
GL eqn. over the disc thickness. Since the order pa- 
rameter does not change in the z— direction, both the 
superconducting current and the vector potential have 



no z— component. Then the boundary condition (Q) is 
automatically fulfilled on the upper and lower disc sides. 

Our 3D calculations show that the discs studied exper- 
imentally in Ref. [1] exist in a regime of surface super- 
conductivity, or the giant vortex state i.e., between the 
Hc2 line and the Hcj, line. If the thickness of the discs 
is further reduced then we go over to the mixed regime 
and the giant vortex breaks up into many vortices if the 
radius of the disc is sufhciently large, even for a type I 
sample. 

Therefore, we consider the situation with a fixed value 
of the angular momentum L for the order parameter 
"^{p) — F{r)exp{iL(f>), when both the vector potential 
and the superconducting current are directed along e^. 
Then eqns. (|^||) are reduced to the following form 
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where Q(x < 1) = 1, Q(x > 1) = 0; A = e^A; i?, d are 
the dimensionless disc radius and thickness, respectively; 
the brackets <> means averaging over the disc thickness 

< /(P) >= lZ%fiz,p)dz/d. 

The magnetic field created by the superconducting 
current in the disc decreases in strength away from 
the disc as a magnetic dipole field: H ^ Be- 
cause of this, in our numerical calculations the condi- 
tion for the vector potential is transferred from infin- 
ity to the boundaries of the simulation region as follows: 
A{z,r = Rs) = \HnRs, A{\z\ = = \Hqp, where 

Rs,ds 3> i?, d are the sizes of the simulation region in 
the radial and z direction, respectively. We use typically 
Rs,ds = (5-^ 10)max(R, d), and checked that an increase 
of the simulation region does not change our results sub- 
stantially. The boundary conditions for the order param- 
eter 
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correspond to zero current density at the disc surface 
and a finite value of the first derivative of F at the disc 
center. To solve numerically the system of eqns. (|[]^) 
we apply a finite-difference representation of the GL and 
3D Maxwell eqns. on the space grid pi, Zj. 

Discs in three different regimes will be considered 1) 
type II, 2) type I and 3) multiple type I behavior. When 
we compare the theoretical results with the experimental 
data we have to keep in mind that experimentally the 
magnetization will depend on the filling fraction of the 
Hall bars used as detectors which is not exactly known. 
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Also because of the square geometry of the Hall detec- 
tor whose sides are of the same size as the diameter of 
the largest disc it will underestimate the magnetization 
(the flux expelled) of the smaller discs. These effects will 
result in an unknown scale factor for the magnetization 
of order 1. In order to have a comparison of relative 
magnitudes such as the size of the jumps in magnetiza- 
tion, we scale the theoretical results such that they have 
the same maximum magnetization as observed experi- 
mentally. When determining the magnetization from the 
LGL eqn., the same method as in Ref. |Q| was used; we 
have assumed the Abrikosov parameter /? to be 1.0 as 
done in Ref. . We compare our theoretical results with 
the experimental results on Al discs at 0.4 K of Geim et 
al 1^ and took for the zero temperature coherence length 
^(0)= 250 nm and the penetration length A(0)=70 nm 



as estimated in Ref. | 
are also given in Ref. 



The disc thickness and radius 
and therefore our theory does 
not contain any fitting parameters. 

Fig. 1 shows the magnetization curves for an Al disc 
of thickness d=0.15 ^.m and radius R=0.315 /xm. Large 
solid dots are the experimental data and exhibit a con- 
tinuous superconducting-normal transition; the dotted 
curve is the solution from the LGL eqn. whereas the 
thin solid curve is the numerical solution of the non-linear 
GL eqns. coupled to the three dimensional Maxwell eqn. 
The dotted curve is scaled by 0.158 and the solid curve 
by 0.537. Surprisingly the dotted curve gives a line shape 
in closer agreement with the experiment but its magni- 
tude is clearly too large. There is some improvement in 
the line shape (dashed curve) if we reduce the disc thick- 
ness to 0.07 /xm in which case the radius of the disc was 
changed to 0.31 in order to keep the critical field the 
same (the magnetization was scaled by 1). Therefore, 
we suspect that the effective thickness of the disc which 
is still superconducting is much smaller than the actual 
thickness. 

Fig. 2 shows the magnetization curves for a larger disc 
of thickness 0.15 /im and radius 0.473 /im at 0.4K. The 
same symbol and curve conventions are used as in pre- 
vious figure. The dotted curve is scaled by 0.124 and 
the thin solid curve by 0.581. It is obvious that the dot- 
ted curve is very different in shape and magnitude and 
shows a jump in magnetization at a very different value 
compared to the experimental curve. This clearly demon- 
strates that a LGL eqn. with a homogeneous magnetic 
field distribution over the disc is not appropriate in this 
case. The finite thickness of the disc results in very im- 
portant geometrical corrections to the field profile which 
influences the superconducting state appreciably. Note 
that the magnetic field dependence of the experimental 
curve is well described by the thin solid curve: 1) the 
slope of the magnetization curve, 2) the non-linear be- 
havior near the step in the magnetization, and 3) also 
the magnetic field at which the step in the magnetiza- 
tion takes place is correctly predicted. The experiment 
shows a first order transition from the superconducting 
to normal state at a magnetic field of 70.86 G while our 



calculations still predict a transition to the L = 1 super- 
conducting state which becomes normal at 81.5 G. The 
origin of this small discrepancy is still not clear to us but 
may be due to effects of disorder . 

The magnetization curves for a disc of thickness 0.15 
/im and radius 1.2 /tm, is shown in Fig. 3. The symbol 
and the curve conventions are again the same as before. 
The dotted curve is scaled by 0.062 and the thin solid 
curve by 0.775. Note that the LGL eqn. in this case 
gives the same type of discrepancies as found in the ex- 
periment of Buisson aZ . Firstly they found that the 
magnitude of the jumps in the magnetization as obtained 
from a solution of the LGL eqn., are too large compared 
to the experimental results. Buisson et al argued that 
this was due to an ensemble averaging in their experi- 
ment. The single disc experiment of Geim et al rules 
out this possibility. It is true that the magnitude of the 
jumps in the single disc experiment is much larger than 
in the many disc experiment, but still, for the single disc, 
the jumps are much smaller than those obtained from a 
solution of the LGL eqn (compare the dotted curve with 
the experimental data in Fig. 3). Our thin solid curve 
gives precisely the same magnitude for the jumps in the 
magnetization as in the experiment and also the correct 
magnitude of magnetization and magnetic fleld for most 
of the transitions. Secondly, they found that the po- 
sition of the first jump in the magnetization obtained 
from the LGL eqn. is much below that of the experi- 
mental curve. No proper explanation could be given for 
this. Our dotted curve gives similar discrepancy of ap- 
proximately the same magnitude and the thin solid curve 
settles this dispute. The first jump coincides with that 
of the experimental curve. We find that if we keep the 
radius unchanged and decrease the disc thickness then 
the upper critical field and the number of jumps in the 
magnetization curve remain unchanged. Only the posi- 
tion of the first peak shifts towards lower magnetic field 
which will of course also lead to an increase in the mag- 
netic field spacing at which the jumps in magnetization 
occurs. This obviously is due to the fact that as the disc 
thickness is reduced the magnetic field inside the disc in- 
creases at a faster rate and so the transition to the first 
fluxoid state occurs at lower applied flelds. Near the crit- 
ical field of course the field inside does not depend on the 
thickness and is the same as the external field. This once 
again shows that the de-magnetization factor is crucial 
for these discs and they determine largely the shape of 
the magnetization curves. Another interesting point to 
be noted is that the experimental curve shows a grad- 
ually decreasing interval of magnetic field at which the 
jumps occur. It was shown in Ref. |Q that within the 
LGL theory flux quantization condition does not imply 
that the jumps in the magnetization will occur at reg- 
ular intervals. As long as the order parameter at the 
central region of the disc is not negligible, the interval 
will decrease slowly and hence only for large values of 
L the interval will be the same as that given by the flux 
quantization condition. But our 3D solution show a more 
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drastic decrease in the interval (since the first few jumps 
do not coincide with the LGL theory whereas the last 
jump does). The reason is that at smaller fields the field 
inside the disc changes very slowly. In fact there is even 
a small regime where the field at the center of the disc 
can decrease with increasing applied field (see curves for 
L=l and L—2 in the inset of Fig. 3). But at high fields 
the field inside the disc increases almost in the same way 
as the external field. In the inset of Fig. 3 we have plot- 
ted the magnetic field distribution H{p, z=0,6=0) of the 
system considered in Fig. 3 for ten values of the applied 
field. The corresponding value of the angular momen- 
tum of the equilibrium superconducting state and the 
magnetic field is also indicated. For L ^ there is sub- 
stantial penetration of the magnetic field in the center of 
the disc while a ring like region near the edge of the disc 
remains superconducting. This is the well known surface 
superconductivity or the giant vortex state that arises in 
a finite sample due to the lower energy of the edge states 
as compared to the bulk states 0. As soon as the sys- 
tem becomes ring like two equal and opposite screening 
currents flow in the system (one along the outer radius 
of the ring like region and the other along the inner one). 
Note that for the highest field shown in the inset of Fig. 3 
there is no bulk superconductivity and only a tiny region 
of surface conductivity survives which is responsible for 
the dip in the magnetic field near the edge of the disc. 

All along we have assumed that the system evolves 
along the free energy minimum and obtained quantita- 
tive agreement for the position, magnitude and period- 
icity of the jumps in the magnetization as well as the 
absolute value of the magnitude of magnetization. But 
there is one discrepancy. Note that the theoretical curves 
in Fig. 3 show a critical field which is appreciably smaller 
than found experimentally and the total number of jumps 
in the experiment is 19 compared to 11 in the theory. 
When we enlarge the disc radius from 1.2 /im to 1.57 
fim, the number of jumps in the magnetization curve is 
increased to 19 but the upper critical field is reduced to 
63 G. It is also to be noted that the slope of the experi- 
mental curve decreases with increasing field and tends to 
become parallel to the field axis at higher magnetic fields. 
In finite systems the Bean-Livingston barrier |]l^ at the 
surface can cause the system not to evolve along the free 
energy minimum. It leads to jumps in the magnetization 
at a much larger value of the magnetic field compared 
to the value at which the L—1 state becomes the ground 
state. We have to decrease the width of the disc to an 
unreasonable value of 0.06 fim in order to have the first 
jump at the same position as that seen in the experi- 
ment for the same disc radius. But, as discussed before, 
such a decrease in thickness does not increase the critical 
field and cannot explain the high field discrepancy. It 
is also known that surface defects can destroy the Bean- 
Livingston barrier for increasing fields. The previously 
mentioned disorder effects ||^ may be responsible for the 
increase of the upper critical field but at present no ex- 
plicit calculation exists with which we can compare and 



therefore this is still an open issue. 
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FIGURE CAPTIONS 
FIG. 1. Magnetization versus the external magnetic field 
for a superconducting disc of radius R= 0.315 /im and 
thickness d=0.15 fim at T=0.4 K. Solid dots are the ex- 
perimental data from Ref. ||] . We show the results of the 
GL theory including the 3D Maxwell eqn. (solid curve, 
and dashed curve for R=0.31 /zm and d=0.07 fJ.m), and 
the result for the LGL-theory (dotted curve). 
FIG. 2. The same as Fig. 1 but now for a disc with radius 
R=0.473 fim and thickness d=0.15 /im. 
FIG. 3. The same as Fig. 1 for a disc of radius R=1.2 
fim and thickness d=0.15 fim. In the inset we show the 
field distribution in the plane through the center of the 
disc for different values of the external magnetic field. 
Far away from the center of the disc the magnetic field 
equals the external magnetic field. 
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